Jamming and growth of dynamical heterogeneities versus depth for granular heap flow 
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We report on grain dynamics versus depth for steady-state gravity-driven flow of grains along 
a heap formed between two parallel sidewalls. Near the surface the flow is steady and fast, while 
far below there is no flow whatsoever; therefore, a jamming transition occurs where depth is an 
effective control parameter for setting the distance from the transition. As seen previously, the 
time-averaged velocity profile along the sidewall exhibits a nearly exponential decay vs depth. Using 
speckle-visiblility spectroscopy (SVS), we find that velocity fluctuations grow relative to the average 
on approach to jamming. Introducing an image-based order parameter and the variance in its 
temporal decay at a given depth, we also characterize the spatially heterogeneous nature of the 
dynamics and find that it increases on approach to jamming. In particular, the important time and 
length scales for dynamical heterogeneities are found to grow almost exponentially as a function of 
depth. These dynamical changes occur without noticeable change in structure, and are compared 
to behavior found by experiment and simulation elsewhere in order to test the universality of the 
jamming concept. In particular we find that the size of the heterogeneities scales as the inertia 
number raised to the —1/3 power, in agreement with two recent simulations. 

PACS numbers: 45.70.Mg, 64.70. Pf 



"Jamming" is defined to occur when a disordered sys- 
tem develops a yield stress, or alternatively, as a prac- 
tical matter, when it develops a stress relaxation time 
that exceeds a reasonable experimental time scale [lj. 
Granular materials are disordered systems that readily 
jam, for example as a static heap that maintains its slope 
indefinitely against gravity 2-4]. The general behavior 
of granular systems remains a topic of fundamental in- 
terest in part because of possible similarities with dense 
colloids, supercooled liquids, and soft glassy materials, 
in their approach to jamming Q. For example the 
microscopic dynamics generally exhibit spatiotemporal 
heterogeneities that are expected to increase in size as 
the relaxation times become longer Such hetero- 

geneities have been directly visualized in dense colloidal 
suspensions Illl - ll3l ] and in granular materials subjected 
to shear [lj-[lg . They have also been visualized for large 
grains in a nearly fluidizing upflow of air [l7j . and have 
been found to increase on approach to jamming with a 
dynamic correlation length exponent [l8|, [ll| in agree- 
ment with diverging length scales observed in simulations 

[Ennui. 

To explore and possibly extend the universality of the 
jamming concept, it is crucial to broaden the range not 
just of systems but also of different control parameters 
that can be varied to adjust the proximity to jamming. 
Toward this end, in this paper, we consider the flow in- 
duced by the steady addition of grains to the top of heap 
confined between parallel plates. Whereas earlier with 
this setup we studied behavior by varying the rate of 
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flow [22j, |23j , now we fix the flow rate and study behavior 
as a function of depth below the surface. In this regime 
of continuous flow, the velocity profile along the sidewall 
decays nearly exponentially with depth [22| ; this observa- 
tion has been reproduced and explored in detail [24T4291 ] . 
In effect the grains are strongly driven and unjammed 
near the surface but become progressively jammed as a 
function of depth. To quantify this jamming transition, 
we measure the coarse-grained surface velocity profiles 
with great dynamic range and we also implement two 
novel techniques to characterize grain-scale behavior. In 
particular we use Speckle- Visibility Spectroscopy (SVS) 
[3Cf to measure the time-resolved relative speed of grains, 
and we introduce a video-based dynamic order parame- 
ter and associated four-point susceptibility to measure 
dynamical heterogeneities, all as a systematic function 
of depth. We find that fluctuations and heterogeneities 
increase with depth, and we compare with similar be- 
havior reported for prior approaches to jamming in both 
experiment and simulation. 



I. EXPERIMENTAL SYSTEM 

The flow cell and granular medium are identical to that 
of Refs. [U HH, in which avalanches were studied by 
slow steady addition of d = {0.25 — 0.35} mm diameter 
glass beads at the top of a heap formed between parallel 
static-dissipating Lucite plates. The width of the gap 
between the plates is w = 0.95 cm, and is closed along 
the bottom and one vertical side; the inner area is square 
28.5 x 28.5 cm 2 . The static granular heap has angle of 
repose 9 r = 24° and bulk mass density p = 1.5 g/cm 3 , 
which corresponds to spheres of mass density 2.4 g/cm 3 
arranged with packing fraction 0.64. 

For all experiments, grains are fed onto the top of 
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FIG. 1: (Color online), (a) Sketch of the experimental ap- 
paratus and coordinate system, plus callout of one raw video 
frame with inverted grayscale. The color code is used in sub- 
sequent figures to indicate depth, (b) Characteristic wave 
number k* in x direction, vs depth z, from the power spec- 
trum of video images, (c) Peak value of the pair correlation 
function g(x), vs depth z, from particle-tracking analysis of 
video images. 



the heap at the back of the box at a fast steady rate 
of Q = 2.5 g/s. This is accomplished using a large funnel 
with a fine control valve, and an aluminum bar sand- 
wiched between the plates to break the free-fall of grains 
and deposit them gently onto the heap. Long-duration 
runs are accomplished by collecting the discharge and 
pouring it back into the funnel. Note that the choice of 
flow rate onto the top of the heap is nearly ten times 
greater than the critical rate Q c = 0.36 g/s [22], [23| . 
below which the discharge out the bottom is intermit- 
tent and accomplished by a series of discrete avalanches. 
Therefore the measurements below are all for a regime in 
which the flow is smooth, continuous, and independent of 
time. Over an extended distance, several channel widths 
in from the entrance and exit, the flow field appears par- 
allel to the top free surface oriented at 9 = 26°, two 
degrees higher than the repose angle. Care is taken that 
all measurements are conducted in this fully developed 
region, such as illustrated in Fig. Q] where the example 
image is acquired about 5 cm from the exit. 

Characteristic scales may be roughly estimated by as- 
suming that the average velocity decreases to zero over 
a depth set by the channel width. This gives a charac- 
teristic flow speed of Q/{pw 2 ) = 2 cm/s, which is nearly 
equal to the projection of grain-size free-fall speed along 
the surface, ^fgdsm9 r . The characteristic strain rate is 
given by speed divided by channel width as 7 = 2 s _1 ; 
therefore, the time 0.5 s to accumulate a strain of one is 
much longer than the grain collision time, which must be 
shorter than \fdfg = 5.5 ms. This means that the grains 
either experience many rapid collisions within a long- 
lived set of neighbors, reminiscent of caging in a glassy 
system, or else they flow subject to enduring contacts. 
Taking pressure as P = pgw, the characteristic "inertia 
number" is I = jd/ y/Pfp — 0.002; this means that iner- 
tial stresses are small compared to confining stresses [|J. 



II. STRUCTURE 

Using a high-speed Phantom video camera, 960 x 480 
pixels and 8 bit gray scale, we begin by acquiring a se- 
quence of still images as shown for example in the callout 
of Fig. [TJ For this, the camera is mounted about 30 cm 
away from the side of the sample, oriented at the same 
angle as the heap, and fitted with a macro lens. The 
magnification is 43.6 pixels/mm, or about 13 pixels per 
grain diameter. Individual grains along the wall are thus 
readily visible, and the intensity profile across grains is a 
series of smooth peaks and valleys (not step functions); 
therefore, their coordinates may be found by standard 
particle-tracking methods. By counting the number of 
particles in a narrow 22-pixel slice for different depths, 
we find that the number density rises from zero to a con- 
stant, presumably random-close packing, by a depth of 
only a couple grains. Two other measures of instanta- 
neous packing structure exhibit similar behavior. The 
first, Fig. Hlb), is for the location k* of the peak in the 
power spectra for image slices at different fixed depths. 
This rises to a constant exactly as for the number density 
profile. The second, in Fig. [He), is f° r the value g* of 
the first peak in the pair correlation function, computed 
from the grain coordinates found in 22-pixel image slices 
at different fixed depths. This, too, rises to a constant 
but somewhat more gradually than k* and number den- 
sity data. By all three measures, the influence of the 
top free surface on the packing structure of the flowing 
grains is limited to a shallow layer of only a few grains, 
much less than the width of the channel. This result 
agrees with numerical simulations for flow on an inclined 
plane [3J Hlj]. Exactly as for the jamming and glass 
transitions, the dramatic changes next discussed for dy- 
namics vs depth into the heap are accomplished without 
significant change in structure. 



III. VELOCITY PROFILES 

The same video set-up and particle-tracking methods 
as above can be used to measure the average speed v x (z) 
of flow in the x-direction along the wall as a function of 
depth z. For this the images are sliced into regions 10- 
pixels wide, and the average velocity at a given depth is 
found by the grain displacement in the x direction be- 
tween successive frames, multiplied by frame rate, aver- 
aged over all particles in the slice, and averaged for 1000 
frames. To measure faster speeds near the free surface as 
well as slower speeds down deep, this is done for a wide 
range of frame rates: {3900, 2048, 1024, 512, 256, 128, 64, 
32, 16, 8, 4, 2, 1, 0.5, 0.25} per second. The results for all 
frame rates are averaged together with weighting given 
by the respective statistical uncertainties. The weight of 
a contribution is highest when the average displacement 
is about 1/10 grain per frame. A similar procedure is 
also employed to measure the average speed v x (y) of the 
flow at the free surface as a function of the distance y 
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FIG. 2: (Color online), (a) Mean speed v x in the flow direc- 
tion, from video data, and fluctuation speed Sv, from speckle- 
visibility spectroscopy data, vs depth z along the sidewall. 
(b) Mean and fluctuation speeds vs distance y along the top 
free surface at z = 0. The depth is redundantly indicated 
using the color code shown in Fig. [T] The gray lines in (a) 
represents exponentials and the and gray curve in (b) is a 
parabola. 



away from the front plate. For this we use 1000 frames, 
576 x 576 pixels, taken at a rate of 4000 per second. 

The mean velocities along the wall and across the top 
free surface are displayed in Fig. [3] The flow profile along 
the wall, Fig. Ufa), begins at 10 cm/s and decays by a 
factor of nearly 10 6 over a distance of about 1 cm. The 
shape appears to cross over between two exponentials at 
about 0.4 cm. Exponential decay has been reported for 
both heaps [H, Hj-Hl] and rotating drums 133: s imilar 
profiles have been observed in Couette cells |34l436j . The 
flow profile across the top surface, Fig. [2fb), also begins 
at 10 cm/s but increases much less dramatically to a 
maximum of 15 cm/s at the center; the shape is nearly 
parabolic. Assuming such a parabolic increase away from 
the wall at all depths, we integrate the granular flux over 
the entire cross section and arrive at a mass discharge 
rate within ten percent of the known value Q = 2.5 g/s. 
Based on this check we suppose that the behavior along 
the sidewall is at least somewhat representative of that 
in the interior. The influence of the side wall on bulk 
behavior was studied in Ref. [27| . 



IV. VELOCITY FLUCTUATIONS 

Next we measure the size Sv of microscopic veloc- 
ity fluctuations by use of Speckle- Visibility Spectroscopy 
(SVS) [30]. This is a time- resolved dynamic light scatter- 
ing method, which we have previously used to measure 
dv vs time for intermittent avalanches on the top of a 
heap 23]. One advantage of SVS over particle track- 
ing for measurement of fluctuations is that the incident 
photons provide an ensemble average over an interior 
region several grains from the boundary. Another ad- 
vantage is that it probes wavelength-scale motion much 
smaller than the grain size; this is particularly helpful 
on approach to jamming where the persistence length for 
ballistic motion becomes so small that particle tracking 
would return an average speed contaminated by multiple 
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FIG. 3: (Color online). Fluctuation speed Sv vs time for 
several depths, z, as labeled. The curves are decorated with 
a solid triangle at the time r/v x needed to advance with the 
average flow by one grain radius. 
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FIG. 4: (Color online). Fluctuation speed Sv vs (a) average 
flow speed v x and vs (b) strain rate times grain diameter. 
Each point corresponds to a different depth, as indicated by 
the same color code as in the previous figures. In each plot the 
solid line represents y — x; in (a) the dashed line represents 
the best fit to a square-root power law. 



direction changes due to neighbor collisions. Here we use 
the same setup as before [23|]: illumination by Nd-YAG 
laser with wavelength A = 532 nm and maximum power 
4 W; detection by 8-bit linescan camera with 2 x 1024 pix- 
els and 58 kHz frame rate. The main difference is that 
the incident light is spread out with a cylindrical lens 
to form a line about 5 cm long and 1 mm wide at vari- 
able depth parallel to the top free surface of the heap. 
The depth-resolution is further smeared beyond 1 mm 
because incident photons migrate within the sample, pri- 
marily over a region set by several grain diameters. The 
camera is placed about 10 cm away, so that the speckle 
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size is closely matched to the pixel size. As the grains 
move relative to one another, the speckle pattern changes 
and hence each image blurs out to an extent that depends 
on the duration T of the exposure. This effect is quan- 
tified by the variance V(T) in detected intensity levels, 
which we compute for a sequence of synthetic exposure 
durations by addition of a given number of successive 
frames of video data. Here, where the laser light is dif- 
fusely backscattered from ballistically moving particles, 
the spectrum is Lorentzian with some linewidth T. The 
theory of SVS [30] gives the variance ratio for exposures 
of 2T and T duration as 

V{2T) _ er Ax - 1 + Ax 
V(T) ~ 4{e- 2x - 1 + 2x) ( ' 

where x = TT ps (4ttSv/X)T. The quantity Sv repre- 
sents an ensemble average over the scattering volume of 
the characteristic speed for the relative motion of grains 
separated by one photon transport mean free path, I*. 
In general this includes contributions from randomly- 
directed fluctuations Sv ran d as well as from shear. For 
pure shear deformation, with no fluctuations, at abso- 
lute rate I7I, the time scale ts in Eq. (7) of Ref. (37j 
gives the expectation for Sv as |7|i*/V5; this is the same 
relationship not ed p reviously for Diffusing- Wave Spec- 
troscopWDWS) HaHH. Since I* is a few grain diameters 
[22I l28l I381 - I40T ] . the relative speed expectation reduces to 
approximately |7|d for pure shear. If the grains addi- 
tionally experience random ballistic velocity fluctuations 
around the mean shear flow, then the relative speed de- 
duced from Eq. (JXJ) will be commensurately larger than 
this lower bound as Sv = Sv ran d + |7 Id- 
Example data for the instantaneous characteristic rel- 
ative speed Sv vs time are displayed in Fig. [3] for five 
different depths. The relative speeds decrease with in- 
creasing depth, and appear to be constant in time. We 
anticipated, but did not find, temporal fluctuations in 
Sv that increase with depth. Such fluctuations could be 
present, but remain unobserved due to the spatial aver- 
aging in proportion to illumination area and the resulting 
photon sampling volume. In either case, we average each 
such time sequence to arrive at a single value for the 
mean relative speed at a given depth. The results are 
plotted in Figs. [2] along with the profiles of mean speed. 
At the top free surface, Figs.EJb), the value of Sv is much 
less than the mean speed. Thus the instantaneous grain 
velocities point primarily in the flow direction, with only 
slight fluctuations around the average. The same is true 
along the wall but only near the top, Figs. [2ja). In par- 
ticular Sv decays with depth less rapidly than v x , so that 
below 0.6 cm the relative grain speed becomes greater 
than the mean flow speed. Below this depth the grains 
don't merely flow but must jostle more and more with 
their neighbors. This can be thought as similar to the 
onset of caging effects on approach to jamming. 

To examine the relationship between fluctuation and 
flow, we display Sv vs v x on logarithmic axes in Fig. EJa) . 
This is a parametric plot, in that each point represents 



a single depth. For the upper two decades in Sv, and 
the corresponding upper four decades in v x , a square- 
root relationship Sv ~ ^JvZ holds quite well as shown in 
Fig-EIa) by the dashed line of slope 1/2. Since Sv and v x 
are not proportionate, there is no quasistatic limit where 
the nature of the dynamics merely slows down without 
alteration of character. Similar behavior was measured 
for continuous flow in a hopper [39( and also at the free 
surface of a heap [23[ , for both of which the flow and fluc- 
tuation speeds were decreased by a reduction in discharge 
rate. The same behavior, Sv ~ y/vZ, is also implied by 
the simulation result Sv 2 oc 7 [41j, since here the veloc- 
ity profiles are nearly exponential; this is seen directly 
in Fig. HJb) where Sv data are plotted vs strain rate. 
In all cases the relative rise of Sv over the flow speed for 
greater depth signals increased jostling and dissipation at 
decreased driving rate. For the hopper this causes clog- 
ging; at the top of a driven heap this causes a transition 
from steady flow to intermittent avalanches; here along 
the sidewall for steady heap flow this causes jamming at 
great enough depth. 

The separate contributions of fluctuation and shear to 
the SVS signal may be judged by the plot of Sv results 
vs |7|d in Fig. Q|b). Since the flow profile is roughly ex- 
ponential, this plot appears similar to Sv vs v x in part 
(a) of this figure. However now we may compare directly 
with the prediction Sv ~ \ j\d for pure shear. Indeed for 
the rapid flow near the surface, the data roughly follow 
this form. But deeper below the surface, the Sv data rise 
above \ j\d. This is an alternative way of demonstrating 
that the grains jostle around the mean flow to an extent 
that is small near the surface but that increases system- 
atically vs depth on approach to jamming. 



V. DYNAMICAL HETEROGENEITIES 

The transition to jamming in glasses, colloids, and 
granular media appears to be universally accompanied 
by the growth of dynamical heterogeneities. Grain-scale 
motion is not possible unless there is cooperation with 
neighbors, for example in the form of string- like swirls 
that come and go; the closer to jamming, the more in- 
frequent the swirls and the larger the number of grains 
involved. A well-developed approach for characteriz- 
ing such fluctuating correlated dynamics begins with an 
ensemble-averaged overlap order parameter Q t (r) vs de- 
lay time r and its average Q(t) = (Qt(i~)) over all start 
times t; see the review in Ref. [42||. For a finite system of 
N grains, spatially heterogeneous dynamics manifest as 
fluctuations in the instantaneous number of fast-moving 
regions and hence as variability in the decay of Qt(r) for 
different start times. This is captured by the the normal- 
ized variance X4( T ) = N[(Q 2 (t)) — (Qt(i~)} 2 ], which is a 
four-point susceptibility that depends on the time inter- 
val t and a spatial length usually comparable to the grain 
size. Note that this definition renders Xi( T ) independent 
of N since the variance for counting statistics decreases as 
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FIG. 5: (Color online), (a) The average overlap order param- 
eter Q, and (b) the corresponding susceptibility \i as func- 
tions of lag time t, for several depths as indicated by color and 
legend. Dashed curves in (a) represent fits to a compressed 
exponential, Q(r) = exp[— (t/to)^], where both to and f3 are 
adjusted. 
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FIG. 6: (Color online), (a) Exponent /3, and (b) decay time 
scale ro, for the compressed exponential fits in Fig. [5ja), ver- 
sus depth z. The time r/v x to flow one bead radius, the 
time r* at which \± peaks, and 1/ |-y j , are also included in 
(b) as labeled, (c) The number n* of fast-moving grains in a 
dynamical heterogeneity versus depth. 



1/N. An explicit counting statistics argument for the av- 
erage and variance of the number of fast-moving regions 
[l8l | then gives the number n* of grains in a dynamical 
heterogeneity as 



Xi 



(Qi - QoXQi - Q*) 



(2) 



where Xi* and Q* are the respective values where Xi( T ) 
is maximum, and where Q\ and Qq are the average con- 
tributions at this time coming from the slow and fast re- 
gions respectively. On approach to jamming, the values 
of X4* an d n* increase because greater cooperativity is re- 
quired for rearrangements [13, EH • If the dynamical het- 
erogeneities are one-dimensional string-like objects, as is 
the usual case, then the growth of n* can be thought of as 
proportional to a growing dynamical correlation length. 

To investigate the possibility of the growth of fluctuat- 
ing dynamics as a function of depth along the sidewall for 
the steady heap flow of previous sections, we employ the 
above formalism using a novel image-based overlap order 
parameter. Recall that the raw image slices obtained ear- 
lier for computation of g{x) and v x have grayscale vari- 
ations indicative of grain positions. It is not necessary 
to identify and track the particle coordinates. Instead, 
much more simply, a suitable overlap order parameter 
may be defined by forming the temporal autocorrelation 
of a sequence of image slices, subtracting the baseline, 
and dividing by the intercept: 



Qt(r) 



{him-m)) 2 



(3) 



where the average is over all pixel indices i in the image 
slice. By construction this decays from one to zero over 
a timescale r/v x set by grain radius divided by average 
flow speed. More generally the length scale probed by 
the "order parameter" Qt( r ) is given by the correlation 
length of intensity fluctuations in the grayscale images, 
which here is about one grain radius. The actual de- 
cay vs t for a given start time t varies according to the 
number fluctuations and the size of the heterogeneities, 
as reflected in the normalized variance Xi( T )- Here the 
image slices are 22 pixels thick and have about N = 150 
grains. 

The experimental results for Q(t) — (Qt(r)) are dis- 
played on linear-logarithmic axes in Fig. [5fa) for a va- 
riety of depths z. In all cases Q(t) decays monotoni- 
cally from one but then exhibits attenuated oscillations 
(not displayed) around zero, similar to the behavior of 
the pair correlation function. Nearer the top free sur- 
face for faster flow, the decay of Q{t) is more rapid; 
deeper down the decay becomes progressively slower, as 
expected. Except for the final oscillations around zero, 
the decay shape may be accurately fit to a compressed ex- 
ponential, Q(t) = exp[— (r/r )' 3 ], as shown by the dashed 
curves. Fitting results for the exponent /3 and the decay 
time To are displayed vs depth in Figs. [HIa,b). Near the 
top surface, the exponent is slightly less than (3 = 2; 
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such a value would be expected for purely convective 
flow, as verified by the spatial intensity autocorrelation 
of individual video image frames. The exponent then de- 
creases with depth and saturates at a value near f3 = 1; 
such a value would be expected if there were considerable 
jostling in addition to flow, as when diffusion dominates 
convection. The fitted decay times To increase monotoni- 
cally with depth by many orders of magnitude. As shown 
in Fig. [6f b) they closely track the time scales r/v x and 
I/7, but are smaller by respective factors of about 2 and 
20. 

The corresponding results for Xi{ T ) vs r are displayed 
on logarithmic axes in Fig. EIb), underneath Q(t), and 
are identically color-coded according to depth. These 
susceptibilities all begin at zero for short delays, rise to 
a peak \i* at a delay time t* that coincides closely in 
Fig.(6jb) with the fitted decay time To, and then fall back 
to zero. Most significantly, the peaks in Xi{ T ) increase in 
height deeper into the sample where the relaxation time 
is slower. Therefore, the dynamics do indeed become 
more unsteady and spatially heterogeneous on approach 
to jamming far below the top surface. 

The size n* of heterogeneities may now be deduced 
at each depth from the peak value \i* using Eq. ([2]). 
For this we estimate the values of Qo and Qi from av- 
erages of Qi(r*) computed separately for pixels i with 
Qi(t*) < Q* and with Q % (t*) > Q*, respectively. This 
gives Qo = 0.2 ± 0.1 and Q\ — 0.7 ± 0.1, where the un- 
certainties represent the uncorrelated spread of values for 
three different depths; the total average is Q* — 0.4±0.1. 
With only slight adjustment of these averages within 
the prescribed uncertainties, the resulting size of hetero- 
geneities at the top free surface is n* = 1, as expected 
far from jamming. All results for n* are plotted vs depth 
in Fig. Etc) on semi-logarithmic axes. The data increase 
with depth by about two orders of magnitude over the 
range of depths measured, appearing to grow without 
bound. The functional form of this divergence is nearly 
exponential, n* ~ exp[z/,z„] where z n = 0.17 cm, as il- 
lustrated by the solid line. As a word of caution, note 
that the value of n* becomes larger than the width of 
the constant-depth image strips; therefore the true het- 
erogeneity sizes may be larger than our estimates. 

Recently Hatano [43| simulated the rheology and mi- 
croscopic dynamics of a three-dimensional granular ma- 
terial under steady shear. From the instantaneous ve- 
locity field, he finds a correlation length that diverges as 
a power law on approach to jamming at zero strainrate: 
£ ~ |7|~ 23 . For string-like dynamical heterogeneities, 
the value of n* scales as £. Therefore to compare directly 
with Hatano's results we plot n* vs strainrate in Fig. [7] 
on logarithmic axes. The results are reasonably well- 
described by Hatano's power-law at intermediate strain 
rates, as shown by the line with slope —0.23; a somewhat 
better fit over the whole range is obtained with exponent 
—0.36. However the meaning of direct comparison is not 
clear since the conditions are not identical: In simula- 
tion the shear is uniform and the pressure increases with 
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FIG. 7: (Color online). The number n* of grains in a dy- 
namical heterogeneity plotted vs strain rate. Each point cor- 
responds to a different depth, as indicated by the same color 
code as in the previous figures. The solid gray line repre- 
sents a fit to n* oc |7|~ ' 23 with exponent fixed according to 
the simulation result of Hatano [43|]. The solid dashed line 
represent the best power-law fit. 

strainrate as P ~ 7 1 ^ 2 , whereas in experiment the shear 
is not uniform and the pressure increases with depth and 
hence decreasing strainrate. In the next section we dis- 
cuss another comparison method. 

Before proceeding, we note that while the results for n* 
vs depth in Fig. Etc) are clearly monotonic and well- fit to 
an exponential for small z, it appears that one oscillation 
cycle around the exponential fit develops at the deepest 
depths examined. Non-monotonic behavior for the peak 
in xt nas been reported for granular [3, |45[ and col- 
loidal [46| systems at high packing fractions where, in 
both cases, some particles do not rearrange or are highly 
constrained. Here, data are taken only where the grains 
rearrange. And structurally, there is no difference in the 
packing configuration vs depth as evident by eye or in the 
measures in Figs. njb,c). Therefore we suppose that the 
oscillation is an artifact of poor sampling statistics due to 
both the slow flow and the large value of n* . Specifically 
at the greatest depths, the speed is roughly one grain 
radius per 1000 s, the strain rate is roughly 10~ 4 /s, and 
the value of n* is roughly 100. According to these orders 
of magnitude, in a 24 hour run, the grains move only 100 
radii and accumulate a strain of only 10. 



VI. DISCUSSION 

We have measured two quantities that describe the 
microscopic dynamics of granular heap flow and how it 
changes with depth: the relaxation times t* w tq w 
0.5r/v x ~ O.O5/7, an d the number n* of grains in a dy- 
namical heterogeneity. On approach to jamming, all of 
these quantities grow without apparent bound. To in- 
vestigate the universality of this jamming transition in 
comparison with other systems, we now seek the dimen- 
sionless parameters that control these divergences at in- 
creasing depth. In general structural rearrangements can 
arise due to stochastic "thermal" fluctuations and also 
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FIG. 9: (Color online). The dependence on reciprocal iner- 
tia number 1/J = ^/gz/^d) of (a) the dynamical times of 
Fig-EJb), divided by the characteristic time d/^/gz, and (b) 
the heterogeneity domain size n*. The solid line in (b) is a 
power-law fit, n* oc (i//)0 33±o.02 



due to shear. The jamming hypothesis is that physical 
quantities ought to be universal functions of the dimen- 
sionless (effective) temperature T = pSv 2 /P and inertia 
number / = jd/^/P/ p. Here the pressure increases with 
depth as P — pgz, since the maximum depth is compa- 
rable to the plate spacing, so these reduce to 



T = Sv 2 /(gz), 
I = jd/y/gz. 



(4) 
(5) 



The latter implies that the relevant characteristic time 
scale is set by pressure and grain size and density as 
t c = dj ' y/gz-i which decreases with depth. The two candi- 
date control parameters are evaluated using experimental 
results from previous sections, and are plotted vs depth 
in Fig. [5] Both are largest near the surface, and decrease 
with depth. Note that T is much less than one, and fur- 
thermore is a factor of 100 or more times smaller than / 



at all depths. Therefore we conclude that T is irrelevant 
and that the behavior is controlled essentially solely by 
/. This makes sense because the stochastic grain motion 
is fueled, ultimately, by the gravity-driven shear flow. 

The microscopic relaxation times, rendered dimension- 
less by comparison with the characteristic time r c = 
d/yfgz, and the dimensionless size scale n* of the dynam- 
ical heterogeneities, are both plotted vs inverse inertial 
number in Figs. EJa,b). This is identical in spirit to the 
Arrhenius plots constructed in Fig. 4 of Ref. 19] for a 
monolayer of gas-fluidized spheres. For the time scales 
in Fig. HJa), note that the dashed line representing I/7 
scaled by t c is identical to 1/1. Since r/v x , t* , and tq 
are proportional to 1/7, all the dimensionless microscopic 
time scales therefore grow linearly with the inverse iner- 
tial number. This makes sense because I is the ratio of 
inertial to confining stresses, which reduces to a ratio of 
time scales. 

The heterogeneity length scale in Fig. I^b) also ap- 
pears to grows as a power of inverse inertia number, but 
with a systematic oscillation as noted earlier in the dis- 
cussion of the n* vs depth plot in Fig. [BJc). The best 
fit is n* oc (i//)0-33±o.02^ ag snown by the solid line. 
In remarkable agreement, Hatano's simulation result (43j 



p r*j j 1 / 2 gives the inertial number as / oc j/y/P ■ 



therefore his result £ 



7' 



-1/4 



7 



3/4. 



can be expressed as 



£ ~ (l//) 1 / 3 . Furthermore, Staron et al. have recently 
found that the length scale of velocity correlations be- 
haves as £ ~ (l//) 3 for gravity-driven granular flow 
down an incline [4 71 ]. 



VII. CONCLUSION 

In this manuscript we have reported on several sig- 
nificant advances. In terms of experimental technique, 
we made the first application of Speckle- Visibility Spec- 
troscopy to heap flow. Furthermore, perhaps more im- 
portantly, we introduced a new method for the measure- 
ment of dynamical heterogeneities based on video image 
correlations. It obviates the need for particle-tracking, 
and could be used for any system in which the pack- 
ing units may be at least partially resolved in images. 
In terms of concept, we recognized that the depth be- 
low the surface of a freely-flowing granular heap may be 
considered as a qualitatively new experimental control 
parameter for the approach to jamming. This comple- 
ments prior work in which jamming is approached by 
variation of temperature or density. In terms of discov- 
ery, we observed that microscopic velocity fluctuations 
grow relative to the mean hydrodynamic flow profile as 
both slow down toward jamming. We also discovered 
that dynamical heterogeneities exist and have size that 
increases as a function of depth. Both such behaviors 
have been noted before for systems with uniform driv- 
ing, but not previously for heap flow. The velocity fluc- 
tuations and dynamical heterogeneities ought to result 
in distinct but measurable local pressure fluctuations, 



8 



which would be interesting to compare with those seen 
in compressed granular media under shear [H, |4!| . Dy- 
namical heterogeneities and their growth with depth for 
heap flow are thus a crucial rate-dependent microscopic 
feature, incorrectly neglected in continuum granular hy- 
drodynamic approaches. It may be an unrecognized rea- 
son for overwhelming focus on flow down a rough inclined 
plane, which is not a more natural or important situation 
than heaps but which allows the imposition of a no-flow 
boundary condition. Our results for microscopic dynam- 



ics are further important in terms of similarity with other 
systems when plotted vs a dimensionless physical control 
parameter. This significantly expands and solidifies the 
universality of the jamming concept. 
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